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ABSTRACT 

In the framework of chiral soliton models we study the behavior of static nucleon proper- 
ties under rescaling of the parameters describing the effective meson theory. In particular we 
investigate the question of whether the Brown-Rho scaling laws are general features of such 
models. When going beyond the simple Skyrme model we find that restrictive constraints 
need to be imposed on the mesonic parameters in order to maintain these scaling laws. 
Furthermore, in the case when vector mesons are included in the model it turns out that 
the isoscalar form factor no longer scales according to these laws. Finally we note that, in 
addition to the exact scaling laws of the model, one may construct approximate local scaling 
laws, which depend of the particular choice of Lagrangian parameters. 
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1. Introduction 

The properties of nucleons in nuclear matter (at finite density ) are of great interest. 
Recently, these properties have been intensively studied [1], in the framework of a simple 
soliton model to the baryon. Such models (supported by the 1/Nc expansion in QCD) 
describe nucleons in terms of effective chiral Lagrangians of mesons. The central idea of 
these investigations is the assumption that the properties of the nucleon at finite density 
are to be obtained just by using different values of the parameters in the underlying meson 
Lagrangian. The relationship between these parameters and their zero density values can, 
for example, be estimated by using QCD sum rules @-|| Here we will not question the 
central assumption but rather will study whether the scaling laws, which hold in the simple 
Skyrme model, remain valid when more complicated but also more realistic soliton models 
are considered. 

The typical results of this approach are the Brown- Rho scaling laws |2|, which hold 
exactly in the semi-classical treatment of the Skyrme model with no pion mass term. These 
express the invariance of the combinations: 

M M* 

A -- 



B = ~^r~ = ~^' (1) 

where the starred quantities correspond to evaluations at finite density. M, qa and (r 2 )j =0 
correspond to the nucleons' mass, axial coupling constant and isoscalar squared radius. The 
quantity f n , which is present in the meson Lagrangian, is the pion decay constant. 

It is natural to ask whether these equations hold in the simplest extension of the Skyrme 
model in which a pion mass term is added. We will show they do hold, provided that, in 
addition, m n scales like 

«r 2 )/=o = <V)i=o ■ (2) 

Now, it is well known that the simple Skyrme model is unable to provide an adequate 
description of several nucleon properties. Examples are: the non-electromagnetic piece of 
the neutron proton mass difference |7J , the proton matrix element of the axial singlet current 
|| and the high energy behavior of the phase shifts in pion nucleon scattering M . All these 
short-comings can be improved in soliton models which contain explicit vector meson degrees 
of freedom. It is therefore natural to study the scaling behavior of static baryon properties 
in a more realistic vector meson modelQ. 

We will show that in this model the quantity B in ([!]) is no longer invariant under 
rescaling. However the combination A in ([!]) will be invariant provided that 

^ = 4 and ^ = 4|, (3) 
m p m* gU g*f* 

hold, as well as some scaling laws for coupling constants in the terms of the meson Lagrangian 
proportional to the Levi-Civita symbol. In (|3|) g is the vector meson coupling constant and 
mp is the p-meson mass. 

*The density dependence of static properties was previously studied in such a model where the 
variation of the meson parameters was adopted from Nambu-Jona-Lasinio model calculations. 
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The general strategy for obtaining these scaling laws is to first construct a universal 
expression for the soliton mass by absorbing as many parameters of the effective meson theory 
as possible into field and coordinate redefinitions. In this process not all meson fields can be 
reparametrized because they are subject to restrictive boundary conditions. Subsequently 
these redefinitions are employed to study the behavior of other nucleon properties when 
varying these parameters. 

The scaling laws mentioned above are exact consequences of the models for any choice 
of parameters. It is also possible to find approximate (local) scaling laws which hold in the 
vicinity of a given parameter. An illustration is provided for the vector meson model. 

Tests of various scaling laws should be available after the RHIC facility is completed. 

2. The Pseudoscalar Model 

Although the scaling laws for dense matter were motivated |1[ from an effective meson 
Lagrangian [jrj], which includes the QCD trace anomaly, these laws essentially follow from 



the simple Skyrme model \ Y2 , [T3| . In order to present this model it is convenient to consider 
the root, £ of the non-linear realization, £ 2 = U = exp(ir ■ ir/ f n ), of the pseudoscalar fields. 
Using the vector and pseudovector objects constructed from £ 

v^d^-gd^ and P, = d^ + ^d^. (4) 

the Skyrme Lagrangian assumes the simple form 

£ Sk = -f tr (p^) + J-^tr {\p»,p v \ \p»y\) . (5) 

Substituting the hedgehog configuration £(r) = exp(zr • rF{r)/2) yields the classical energy 
functional 

E cl [F] = 2n^ J™ d({( 2 F' 2 + 2sin 2 F + sin 2 F (<2F' 2 + ^f^j J , (6) 

where the dimensionless variable ( = ef n r has been introduced. This functional is the 
building block for discussing the variation of static nucleon properties with the parameters 
f n and e. The chiral angle, F(Q is obtained by minimizing -E c i[-F] subject to the boundary 
conditions F(0) = —n and F(oo) = for topological baryon number one. Within this 
treatment one derives scaling formulas for the classical mass, the isoscalar radius, (r 2 )/ =0 , 



and the axial vector charge, [P~3 



£ cl = 72.9^, (r 2 ) /=0 = ^ and 9a = ^- (7) 



It should be noted that in the simple Skyrme model the isoscalar radius receives its sole 
contribution from the topological baryon density. This current can be interpreted as the 



Uy(l) Noether current from the Wess-Zumino term [13 



wz 



240vr 2 Jm 5 



f tr(p 5 ) , p = p^dx» , (8) 

J Mr, 



where the differential forms notation has been adopted. Furthermore M 5 refers to a manifold 
which has Minkowski space, M4, as boundary. We are assuming the three flavor case in 
writing (H). 
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In leading order of the 1/Nc expansion the nucleon mass, M, is identical to E c \. Hence 
by eliminating the Skyrme parameter, e, from the expressions ((7|) one easily obtains two 
quantities which only contain physical observables 

M (r 2 \ - f 2 

A := — — = 17.2 and B := V /J ~° J " = 0.062 , (9) 
J-kJIa 9 a 



where the data refer to the Skyrme model predictions ([/]). In particular, these combinations 
are universal in the Skyrme model, in the sense that A and B remain unchanged when 
scaling the fundamental parameters of the Lagrangian (Q) [I], [TJJ 

A* = A and B* = B . (10) 

Here the asterisk denotes quantities, which are computed using^ /* ^ f n and e* ^ e. 

It is an easy exercise to verify that A and B are universal quantities independent of 
which higher order (in derivatives) stabilizing term is adopted. The only requirement is the 
restriction to a single stabilizing term. 

As the Skyrme model represents the simplest version of an effective meson theory which 
supports soliton solutions, it is certainly incomplete. Hence the natural question arises 
whether the scaling laws (|10"D hold in more realistic models. The first extension, which 
comes into mind in this context, is the inclusion of a pion mass term. This neither changes 
the analytical form of (r 2 )j =0 nor of g^] however it contributes to the soliton mass: 

POO 

(^U/e)(mJefA 2 / d^\l - cosF) . (11) 
J o 

Hence a universal expression for the nucleon mass can only be obtained by demanding that 
m^/ef^ does not scale. This implies 

<V>i=o = ™l(r 2 ) I=0 ■ (12) 

We thus see that the scaling laws flTDD do not simply follow from a chirally symmetric 
theory; rather their implementation yields relations among the scaling behaviors of various 
hadron observables. Note that this relation goes beyond dimensional analysis since the model 
contains two dimensional parameters: f n and m n . From Nambu-Jona-Lasinio model studies 
on the density dependence of m n |l6j one might conclude that (r 2 )* I=0 does not vary until 



three times nucleon matter density is reached. In any event, the role of the pion mass is 
special as it arises from the explicit breaking of chiral symmetry. In the following sections 
we will therefore discuss such constraints for coefficients of chirally symmetric expressions in 
more realistic models which include vector mesons. 



3. The Vector Meson Model 

The chirally symmetric action for a realistic vector meson model has been worked out 
previously |I7|, [L8| 

A = J d 4 xC nan + T an + T wz (13) 

^Although we do not specifically have in mind to restrict the discussion to the density behavior of the 
nucleon properties we follow the conventions of dense matter literature and denote quantities which are 
associated with modified meson parameters by an asterisk. 
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c 



nan 

r 

- 1 - an 



tr 



-f Vtf? + ^(u + U^-2)- l -F, u {p)F^{p) + m\R,R» 



tr 



7i + 



V - -</ 7a F(p) [p^R - Rp] - g 2 [72 + 2 73 ] ) • 



For convenience the differential forms notation has again been used to simplify the presen- 
tation of the e pvpa terms. The vector mesons ui and p are contained in R p = p p — (i/2g)v fl 
with Pn = (o> M + p"r )/2. The additional parameters g and 7, are related to the decay widths 
of the vector mesons. Using a slightly different notation they are found to be |18| 



7 2V2 _ 7s 

g « 5.6 , /i = — 71 ss 0.4 , g vv<t> = g l2 ~ 1.9 , « = — « 1 

3 72 



(14) 



The value for k was obtained from a fit to static nucleon properties. 

In addition to the hedgehog ansatz for the pseudoscalar fields we have two more radial 
functions parametrizing the vector meson fields 



W (r) _ G(r) 

, Pi 



gm v 



gr 



(15) 



where my = V^gf^. Introducing a scaled radial coordinate ( = m v r yields the classical 
energy functional [IF 



cl 



my 



jf dc {l ( F ' 2(2 + 2sin2F ) + \ Ae (1 ~ C0SF) ~ ^' 2 + (2 



G 



+G + ( G + 2 Y + Vpi 1 + G ~ cosF ) + 2c/7iF'cjsin 2 F - 4^ 72 G'a;sinF 
+2g-f 3 F'ujG (G + 2) + 2# ( 72 + 73) F'cu (l - 2(G + l)cosF + cos 2 F) I. (16) 



Here a prime denotes the derivative with respect to ( and pi = mi/my are dimensionless 
mass parameters. In order to obtain a universal mass functional we demand that none of 
the coefficients in eq ([II]) vary when f n and/or g are scaled. This implies 



ml 



m. 



m. 



9* ft gfn g*f* gfn 



g *h* = gh , g VV(j> = g V v<t> and k* 



K . 



(17) 



The second equation states that for a universal mass functional the KSRF |H| relation has 
to be taken invariant. For g* 7^ g the frequently adopted scaling relation |T], |2T], 
m*/m p = f*/ f n is apparently violated. 

Static properties are computed as the appropriate matrix elements of the symmetry 
currents. In order to obtain these currents one first introduces external left and right gauge 
fields (Bl and B^) such that the action (13|) is invariant under local chiral transformations 
(up to the anomaly). This introduces an additional term, which does not contribute to the 



classical energy functional, in the anomalous sector [l8j 



d x J tr (£F(B R )£ + ?F(B L )£) (Rp - pR) 



18) 
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This part of the gauged action together with the 7 2 -term contributes to the decay width 

g VVlj) /2g 2 \ 



23 . The currents 



of the process uo — > 7r°7. This determines \2d\ — gyy^j2g l \ fa 0.038 
are extracted from the expression linear in these gauge fields. Substituting the static field 
configurations into the so-obtained currents and taking matrix elements provides the static 



nucleon properties. For the axial charge one finds |p3 

47T 

9g 2 JO 

sinFcosF + 2/i^sinF (1 
cosF (wC-w'tt 



9a 



- 2 j Q d C (2(a 1 + a 2 ) 



(19) 



-912 



cosF) + g (27! + 37 2 ) ujF' sinFcosF 

2 



0-2 



G 

G - cosF)) - uo'sitfF 
+g (72 + 273) F'usmF(l +G- cosF) + Ag 2 d x (a/sin 2 F + uuF'sm.2F 

C 2 F' + g (2 7l + 3 72 ) ^sin 2 F - g^G(G + 2) 
+g (72 + 273) u(l + G- cosF) 2 + 8g 2 d 1 usm 2 F . 



Although the term proportional to d\ is a total derivative and hence does not contribute 
to gA, it is clear that demanding g* 2 d\ = g 2 d\ in addition to the scaling laws (O) enables 
the higher moments of the axial current to scale universally. Note that the first Brown- 
Rho scaling condition in (|1(]) A* = A is satisfied. Starting from the expression for the 
pion-nucleon coupling constant [p3]l 

9-kNN = ^rMUw 2 / dr rhmF{r) (20) 
9 Jo 

it is straightforward to verify that g* NN f*/M*g A = gnNNfw/MgA when (|17|) is imposed, i.e. 
the Goldberger Treiman relation is scale independent. From (^) it can be seen that the 
Wess-Zumino term (|8|) does not contribute to g^- This is in contrast to the isoscalar radius 



23 



(r 2 )i- 



S7T 



m v N c Jo 



d(( 2 R 



1=0 



(21) 



R 



1=0 



C 2 co + 



c 



Att 2 



7i + ^72 



F'sin 2 F + — [F'G(G + 2) - 2G"sinF] 



29 



72 + 273 
29 



F'(l + G- cosFf 



d 



+ Ad x — [sinF(l + G — cosF)] 



where the term which involves Nc stems from Twz- This term scales differently from all the 
others and causes (r 2 )i =0 not to scale universally. If the Wess-Zumino term had yielded the 
sole contribution the isoscalar radius would scale like l/m v ~ l/g*f*. As a consequence the 
quantity B, which is defined in eq @, would be scale independent and one would recover the 
simple Skyrme model result. Also the relation (0) would remain valid. On the other hand, 
if the other terms, which represent the vector contributions, were the only ones present the 
scale independent quantity would instead be 



B :-- 



(r 2 )i- 



9 2 a 



In figure 1 we display the numerical comparison of the scale dependences of B and B. 
this study the variations of the meson parameters are described by the ansdtze 



(22) 



For 



/; = /.a 



and 
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g(l - ex) 



(23) 
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According to eq (\L7\) this implies m* ps m p (l — (1 + c)x) for small x. With c = —0.34 we 
obtain m*/m p ps 0.78 and g\l 9a ~ 0.80 at x = xe = 0.35 as suggested for nuclear density by 
QCD sum rules || || and studies of the neutron beta-decay in heavy nuclei p4| , respectively. 
Note that c 7^ is necessary for qa to vary. The numerical calculation yields M* /M m 0.60 



0.00 o.os 0.10 0.1 s 0.20 



Figure 1: Comparison of the two scaling behaviors for the isoscalar radius in the vector 
meson model. These quantities are normalized to their unsealed values. The variation of 
the soliton properties corresponds to c = —0.34 in eq 



at x = xe, which is slightly smaller than the QCD sum rules estimate of 0.67 ± 0.05 [§, [§]. 

The scaling behavior shown in figure 1 indicates that B is closer to an invariant than B\ 
the variation of the former is only about 5% compared to a 17% decrease of B at xe- The 
physical interpretation of this result is that the vector mesons provide the major contribution 
to the isoscalar form factor. This in turn alters the scaling law obtained in the simple Skyrme 
model which does not contain any vector meson degrees of freedom. On the other hand the 
isovector radius (r 2 )j = i has no direct contribution from the Wess-Zumino term. Using 
formulae (3.5a) and (B2) of ref p3| one verifies that (r 2 ) 1=1 f% / ' g a does not change when (|PTD 



is imposed. This behavior is also obtained in the pure pseudoscalar Skyrme model. 

For completeness we should note that a similar behavior, i.e. the non-existence of a 
universal scaling law is also observed for the isoscalar magnetic moment. Hence the above 
discussion applies to the whole isoscalar current. 

These studies raise the question whether it is possible to construct a vector meson model 
without the Wess-Zumino term. Noting that this term is crucial not only from a conceptual 
point of view (chiral anomaly) but also for the proper normalization of the nucleon charge this 
question can immediately be answered in the negative. It turns out that the normalization of 
the isoscalar charge does not completely fix the scaling law for the associated radius because 
the spatial integral over those parts of the isoscalar density which do not arise from the 
Wess-Zumino term vanishes identically P3 . 
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4- Local Scaling Laws 



The scaling laws (17) for the mesonic parameters, which were imposed to obtain a uni- 
versal mass functional, are clearly very restrictive and one wonders whether or not other 
relations can be found. While the scaling behavior of the parameters entering £ n an acquires 
some justification from QCD sum rules there is no a priori analysis which prevents one from 
choosing scaling laws different from (|TT| ) for the parameters describing r an . In turn a dif- 
ferent choice might yield scaling laws for the soliton properties which significantly deviate 
from those obtained in the simple Skyrme model. In contrast to the general relations given 
in eq ([T7|) such alternative scaling laws will depend on the particular values the mesonic 
parameters take on at zero scaling. We therefore refer to these relations as local scaling laws. 
As this subject is potentially vast we will consider only one example: 



h* = h(l - cx) and g* vv<t) = g V vA l ~ cx f 



(24) 



while /* and g* are taken as given in eq (23). The unsealed parameters are given in eq ([b 
The associated Brown-Rho scaling as well as a modified Brown-Rho scaling, A* = M* / f*g% 
are shown in figure 2. Although the ^-dependences of both quantities are very moderate, 



M7f; g ; 





0.10 

i-fX 



0.10 

i-fX 



Figure 2: Comparison of the two local scaling laws in the vector meson model. 

A* apparently varies less than A*, cf. eq @. Certainly a fine-tuning of either the scaling 
rules (^) or the modified Brown-Rho law could be imposed to yield A* invariant. This 
simple example of a local scaling law indicates that other scaling behaviors for the baryon 
properties could be justified by a different choice for the variation of the parameters entering 
r an . Hence the confirmation of specific scaling laws for static nucleon properties from a 
realistic soliton model seems to be difficult as long as the a priori scaling behavior of the 
anomalous terms is as poorly known as at present. When axial vector mesons are included 



25] the situation is most likely even more arbitrary. 
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5. Conclusions 

In the framework of chiral soliton models we have studied the dependence of static nucleon 
properties on rescaling of the parameters describing the effective meson theory. The first step 
in deriving scaling relations between nucleon properties and the mesonic parameters consists 
of constructing a universal energy functional for the soliton. In the case of the simple 
Skyrme model this immediately leads to parameter dependences of the axial charge and 
the isosinglet radius which obey the Brown-Rho scaling laws. When extending the Skyrme 
model by adding the pion mass term and/or higher order stabilizing terms, the scaling laws 
no longer follow automatically. However, it is possible to impose conditions on the additional 
parameters such that the scaling laws are recovered. In the pseudoscalar Skyrme model the 
isoscalar form factor completely arises from the Wess-Zumino term, which has no effect on 
the classical mass functional. Therefore the conditions on the mesonic parameters, which are 
derived from the mass functional, have no direct consequences for the isoscalar form factor. 
The situation changes drastically when explicit vector meson are contained in the soliton 
model. Then the isoscalar form factor receives contributions not only from the anomalous 
Wess-Zumino term but also from the direct coupling of the photon to the u meson and its 
source terms which are contained in the chirally symmetric e-terms. Since these quantities 
also appear in the mass functional there is no unique scaling law for the isoscalar form factor. 
Numerical calculations indicate that the vector terms actually dominate the isoscalar form 
factor, making necessary a modification of the Brown-Rho scaling law for the isoscalar radius. 
In the context of the vector meson model we have also seen that a universal mass functional 
requires the KSRF relation to be scale invariant rather than just the dimensionless ratio of 
the vector meson mass over the pion decay constant. In addition to this ratio the KSRF 
relation involves the vector meson coupling constant. This coupling constant has to be taken 
scale dependent in order to allow the axial charge to vary. 

These studies were motivated by the derivation of the Brown-Rho scaling for the de- 
scription of dense matter from the simple Skyrme model. In general this approach is based 
on the assumption that the density dependence of nucleon properties can be obtained by 
a suitable scaling of the parameters contained in the mesonic Lagrangian from which the 
soliton is constructed. Unfortunately only little is known about the density dependence of 
the coefficients of the e terms. As these play an important role for stabilizing the soliton it is 
not surprising that by choosing an arbitrary behavior of these parameters almost any scaling 
law for the nucleon mass can be obtained. We have illustrated this for the case that Mj f n gA 
can be made less sensitive to parameter variations than Mj f ny /gA. When applying these 
results to investigate the density dependence of nucleon properties one should also bear in 
mind that at finite density additional terms may be needed in the effective meson theory as, 
for example, suggested by heavy baryon chiral perturbation theory f26f . 
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